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Abstract 

We examine homogeneous but anisotropic cosmologies in scalar-tensor 
gravity theories, including Brans-Dicke gravity. We present a method for 
deriving solutions for any isotropic perfect fluid with a barotropic equation 
of state {p oc p) in a spatially flat (Bianchi type I) cosmology. These 
models approach an isotropic, general relativistic solution as the expansion 
becomes dominated by the barotropic fluid. All models that approach 
general relativity isotropize except for the case of a maximally stiff fluid. 
For stiff fluid or radiation or in vacuum we are able to give solutions for 
arbitrary scalar-tensor theories in a number of anisotropic Bianchi and 
Kantowski-Sachs metrics. We show how this approach can also be used to 
derive solutions from the low-energy string effective action. We discuss the 
nature, and possibly avoidance of, the initial singularity where both shear 
and non-Einstein behavior is important. 
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1 Introduction 



Scalar-tensor theories of gravity [0, ^ allow the gravitational coupling to vary, 
becoming a dynamical field rather than a fixed constant. This occurs in a range of 
fundamental theories that seek to incorporate gravity with the other interactions. 
In Kaluza-Klein models this arises from the variation of the size of the internal 
dimensions. In string theories the dilaton is a scalar field that is necessary for 
the consistent description of the motion of a string in a curved spacetime [Q. 
In a cosmological context such theories allow one to seek dynamical answers for 
questions such as why the Planck mass is so much larger than other physical 
scales or why it should appear to be fixed today. 

The general form of the extended gravitational action in scalar-tensor theories 

is 



S 



/ d'^X\/—q 

167r y ^ ^ 



[I) 



167r 

where in general relativity mpj = F{lp) remains a constant. Brans and Dicke's 
model of gravity corresponds to the particular choice of [/ = and -^(v^) = 
!f'^/8uj where u is a. constant parameter. The more general Lagrangian can still 
be re-written in terms of their Brans-Dicke field = F{(f) with the Brans-Dicke 
parameter uj{(p) = F/[2{dF/d(p)^] 



(PR - ^(V0)2 - + 167r£^atter 



(2) 



and it is this form that we shall use here with U{(f)) set to zero to ensure a strictly 
Newtonian weak field limit. The post-Newtonian parameters of general relativity 
are also recovered in the limit that u ^ oo and {(j)/uj^){duj/d(f)) — > The case 
where U{(f)) is non-zero but = is equivalent to higher-order gravity theories 
[^] with Yukawa type corrections to the Newtonian potential. 

Most analytic cosmological solutions until recently were restricted to the case 
of Brans-Dicke gravity or a few other specific choices of Lj{(f)) ||, ^, |10], [Tl|, [12|. 



Here we present techniques for deriving anisotropic cosmological solutions for 
general scalar-tensor theories where u may be an arbitrary function of 0. This 
can lead to radically different evolution of the Brans-Dicke field as well as the 
Brans-Dicke parameter. Throughout we draw heavily on results derived in the 
conformally related Einstein frame [1^, introduced in Sect. @. In Sect. 



we define the quantities we will use to describe the evolution of homogeneous 
spacetimes. 

We give solutions for barotropic fiuids, including dust and a false vacuum 
energy density in Sect. (El), extending to Bianchi type I models the method used 



recently by Barrow and Mimoso flj] in spatially fiat FRW models. This also 
allows us to consider the role of general relativity as a cosmological attractor 
within a large class of uj{(j)) theories. 

In Sect. (H) we extend the method introduced by Barrow to derive solutions for 



Friedmann-Robertson- Walker (FRW) in vacuum or containing radiation [T^ or 
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stiff matter |TH| , to solve for homogeneous spacetimes including spatial curvature 
as well as shear. In particular we give general solutions for arbitrary uj{(j)) theories 
in vacuum, or with radiation or stiff fluid, in Bianchi types I and V, as well as the 
general vacuum or stiff fluid solutions in Bianchi types III and locally rotationally 
symmetric (LRS) type IX and Kantowski-Sachs, by exploiting known solutions 
in general relativity. 



2 Conformal frames 

Our field equations, obtained by varying the action in Eq. (|]) with respect to the 
metric and field (p, are 

ab -'-„a^)D^ A Q^rpab , ( „ac hd ^ „ab cd\ ^4',c4',d 



+ {g'^'g" - g'-'g'') VcV,0 (3) 

(3 + 2cu)n0 = SttT - (^-^V.^^fe (4) 

where the energy-momentum tensor T"'^ = 2(5£matter/^5'ab)- 

In the scalar-tensor gravity theories the principle of equivalence is guaran- 
teed by requiring that all matter fields are minimally coupled to the metric Qab- 
Henceforth we will refer to this as the Jordan metric. Thus energy-momentum is 
conserved: 

V»T,, = . (5) 

However scalar-tensor theories can be re-written in terms of a theory with a 
fixed gravitational constant with respect to the conformally related "Einstein" 
metric: 

dab = G'Pgab , (6) 

where G is in fact an arbitrarily chosen constant which becomes the gravitational 
constant in the conformal metric. Note that if =const then the two frames 
are identical (allowing for an arbitrary constant rescaling of coordinates) and so 
the scalar-tensor results must be the same as in general relativity whenever this 
occurs. Notice also that for < we must pick a negative G to maintain a 
positive conformal factor. Henceforth we shall assume > 0. 

Instead of appearing in the gravitational Lagrangian, the Brans-Dicke field 
now appears as a scalar field interacting with matter. 



S= I d^x 

1677 
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T^R- o 75 ^ lOTrLmatter 

G 2 0^ 



(7) 



where the covariant derivative Va is taken with respect to the Einstein metric 
and 

r 

p •'—matter /q\ 

i-matter = • l»j 
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Thus, although we recover the famihar Einstein field equations, energy-momentum 
is no longer conserved independently of the Brans-Dicke field 

V'^T,, = -\^f: , (9) 

except when the energy-momentum tensor is traceless — corresponding to vac- 
uum or radiation. In general it is the difficulty of including this interaction 
between matter and the Brans-Dicke field which limits our ability to produce 
analytic solutions. 

Of course the overall energy momentum tensor is conserved (as guaranteed 
by the Ricci identities) as long as we include the Brans-Dicke field as a matter 
field with energy-momentum tensor 

f^' = ir't' - ir'r') ^,c^,d , (10) 



where we define 



_ 3 + M0) 



Here we take 3 + 2u; > 0. [| 

This allows us to deal most easily with other massless fields in the matter 
Lagrangian, exploiting known solutions in general relativity to produce scalar- 
tensor counterparts |jT6|. Short-wavelength modes of a field, oc exp(zfc(ix") 



where kji"' = 0, act like radiation with a traceless energy-momentum tensor. 
Long- wavelength modes of a massless scalar field, v^(t), act like a stiff ffuid with 
density equal to pressure equal to 0^/2 fl^. Although this minimally coupled 
scalar-ffeld in the Jordan frame interacts with the Brans-Dicke field in the Einstein 
frame, they combine to give the same dynamical effect as a that of single stiff 
ffuid. We will also consider massless fields which occur in the low energy effective 
action of string theory. 

We will present in the Sect. (^ results for scalar-tensor gravity in a number of 
anisotropic cosmologies in vacuum, with a stiff ffuid and, in some cases, radiation. 
Before that in Sect. (^ we will give solutions for other barotropic ffuids but 
restricted to a Bianchi type I metric. 



3 Homogeneous spacetimes 

Spatially homogeneous spacetimes admit a group of isometrics acting transitively 



on their spacelike hypersurfaces |]T8[. It is possible to write the metrics of these 
models as 

ds^ = -dt^ + Jaf3{t)W'^aW'^bdx''dx^ , (12) 



-'^If this is not so then the scalar field in the Einstein frame has a negative energy density 
and the Minkowski vacuum, for instance, will be unstable. This seems to be a strong physical 
argument for rejecting such models and henceforth we shall assume 3 + 2a; > 0. 
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where the W^a (dual to W^a) are the invariant vector fields of the reciprocal 
group, [Wa"', Wis""] = apW^"" . The structure constants of the isometry group 
C"^Q,/3, satisfy the Jacobi identities C'j'^^C^]^ = 0. 

These models fall into two classes. Those of Bianchi type and Kantowski-Sachs 
models. The former have a three dimensional group acting in a simple-transitive 
way on the spatial hypersurfaces. They separate into equivalence classes accord- 
ing to Bianchi's classification of the three dimensional Lie groups pO|, |18|. The 



latter has a four dimensional isometry group, but no three simple-transitive sub- 
group acting on the three space. Instead it acts on a two-dimensional subspace. 
This subspace has constant (positive) curvature. 

Although in all spatially homogeneous models the unit normal to the spatial 
hypersurfaces is a geodesic vector field (with tat"" = — 1) invariant under the 
group |T8|, the matter fiow might be tilted relative to this direction. For the sake 
of simplicity, in what follows we shall only be concerned with those models in 
which the velocity of matter is parallel to the unit normal. It follows that T"^ has 
a timelike eigenvector and, thus, both rotation and acceleration are vanishing. 

The homogeneous three dimensional spatial hypersurfaces thus have a metric 
habi orthogonal to the unit vector field, such that the full four-dimensional metric 
in the Jordan frame 

Qab = hab — tah ■ (13) 

A perfect co-moving fiuid, with density p, isotropic pressure p, and (possibly) 
anisotropic stress vr^ then has an energy-momentum tensor 

T^^ptj' + phl + nl. (14) 

We will only consider isotropic matter in what follows with vr^ = (although 
one should bear in mind that an anisotropic expansion can induce anisotropic 
pressures in some fiuids [0). The extrinsic curvature of the hypersurfaces can 
be defined as ^ 

ta;b = -Ohab + (Tab , (15) 

where the expansion 6 = t""^ and the shear = We will also find it 

convenient to define a volume scale factor V such that 

Note that if we redefine the contribution of the derivatives of the Brans-Dicke 
field to the right-hand side of the scalar-tensor equations of motion in Eq. (^ to 
be an effective energy- momentum tensor Tab, so that 

Rab-lgabR)(p = 8n(Tab + fab) , (17) 



2^ 

then we find that the non-minimal coupling of the scalar-field to the spacetime 
curvature creates an anisotropic stress which is proportional to the shear: 

»if ^ 3^^.. . (18) 
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Thus though the Brans-Dicke field will not produce shear where none exists, it 
does exert an anisotropic pressure in the presence of shear even when the other 
matter is isotropic. On the other hand in the Einstein frame the definition in 
Eq. ( ]T0| ) of the energy-momentum tensor Tab, associated with the Brans-Dicke 
field can exert no anisotropic stress for a homogeneous field. 

It is straightforward to calculate conformally transformed quantities in the 
Einstein metric where we find 

2 ;^2_ 



V = (G0)3/V , dt^ = {G(j))dt- , a 



(G0) 



^ ' ^ ■ ^^^^ 

The metric field equations (Q) then yield the familiar constraint equation 

= 247rG (p + p) + Scr^ -^^^^R, (20) 

where ^^'R is the curvature scalar of the hypersurfaces of homogeneity, and the 
Raychaudhuri equation 

dO 1 ~ 

-^+-6^ = -A7rG{3p + 3p + p + p)-2a^ . (21) 
dt 6 

The total energy-momentum in this frame is conserved: 

^ip + p) + e{p + p + p + p) = 0, (22) 
dt 

where we include the density p and pressure p of the Brans-Dicke field. 

But we are also interested in how the Brans-Dicke field and other matter evolve 
separately. The continuity equation for isotropic matter minimally coupled in the 
Jordan frame remains simplest in that frame 

^ + Hp + p)=0, (23) 

which can be integrated for a barotropic fiuid with p = (7 — l)p to give p oc V~^' 
or p oc 

(G0)(3^-^)/2y-7 in the Einstein frame. 
The equation of motion for the Brans-Dicke field is then 



d"^^ -dip 1 / levrG 

This is the one equation of motion in the Einstein frame that still includes an 
explicit dependence on the form of In vacuum or with radiation the right- 

hand side is zero and so the evolution of both ip and the volume scale factor 
V in the Einstein frame is independent of the choice of function uj{4>)- We will 
show in Sect. (^ that these solutions for V{t) also include all the scalar-tensor 
cosmologies containing stiff matter. Analytic solutions for other barotropic fluids 
but restricted to spatially flat FRW models have recently been given by Barrow 
and Mimoso . We will now extend this method to the case of Bianchi type I 
metrics. 
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4 Barotropic matter in Bianchi type I 



We will first consider the case of barotropic fluids (where p = (7 — l)p), such 
as dust (7 = 1) or false vacuum (7 = 0), but excluding for the time being the 
exceptional cases of stiff fluid (7 = 2) or radiation (7 = 4/3) as well as vacuum. 
We will also restrict our analysis to the simplest case of a Bianchi type I metric, 
where the structure constants C^^ all vanish and the spatial curvature is zero. 
The metric can then be written as 

ds"^ = -dt^ + al{t)dx^ + al{t)dy'^ + al{t)dz^ . (25) 

which includes the spatially fiat FRW metric when The expansion 

is given by 

9 - 5| . (26) 

a dt 

where the volume scale factor 

V = = 010203 . (27) 



The '?/'- field equation of motion, Eq. (P^, is driven by the barotropic fluid 
density in the Einstein frame 

3 M(G0)(^^-^)/^ 
^=8^ ' ^^^^ 

where M = SnGpa^'^ /3 is a constant and S'^ = is the volume scale factor in the 
Einstein frame. This same energy density drives the evolution of the three scale 
factors whose individual expansion rates 9i = cii/di, obey the field equations 

^ a 

However the difference between the expansion rates in any two directions is 
not driven by the isotropic fluid. The shear is 

so in the conformally transformed Einstein frame it evolves like the energy density 
of a minimally coupled stiff fluid 

+ 290-^ = . (31) 

dt 

Thus it behaves exactly like a stiff fluid with density 

52 



SttG 327rGa6 



(32) 
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where is a constant, evolving in a flat FRW metric with scale factor a. 

Note that this general relativistic result is not in general true in scalar-tensor 
gravity. Just as a minimally coupled scalar field in the Jordan frame interacts with 
the field ip in the Einstein frame, so shear, which evolves freely in the Einstein 
frame, is coupled to the Brans-Dicke field back in the Jordan frame, and we have 

= 4(G0)2a6 • ^^^^ 

This is a result of the effective anisotropic pressure (in the presence of shear) 
induced by the Brans-Dicke field in the Jordan frame, but absent in the Einstein 
frame. 



We will use the approach developed by Barrow and Mimoso for the scalar- 
tensor solutions in spatially fiat FRW metric, itself an extension of the method 
used by Gurevich, Finkelstein and Ruban [§] for Brans-Dicke gravity. We intro- 
duce the time coordinate ^ defined by 



dt = s3(^-i)(G0)(^-3^)/2 /3 + 2a; ^ ^g^^ 



and the variables^ 



y = V^«^' (35) 

Zi = d% , (36) 

z = \T.^^■ (37) 

i 

The equations of motion, Eqs. ( |2^ and (^), reduce to 

y' = M (4 - 3 7) , (38) 

/ 3 + 2a;3(2-7)M 
= ^ , (39) 



plus the constraint Eq. (pp]), which becomes 

Az^ = y2 + + 4M(G'0)(=^^-^)/2 ~3{2-^) _ ^4Q^ 

Each term on the right-hand side is non-negative and so we can describe the 
expansion in the Einstein frame, 2;^, at any time as being dominated by either 
the term arising from the Brans-Dicke field energy density y"^, the shear S^, or 
the matter energy density 4M(G0)(37-4)/2 53(2-7) _ 

^Note that these variables correspond to those defined by Barrow and Mimoso in terms of 
Jordan frame variables where z = ^(3 + 2aj)/3[(a;/3) + iy/2)] with G = 1 in that paper [p^. 



8 



These equations of motion can be integrated to give 
y = (4-37)M(e-e*) , 



3(2 -7)M 



(41) 
(42) 



where C,o, ^* and o"j are integration constants. Henceforth we will set the constant 
to zero without loss of generality. This merely amounts to a translation of the 
time origin. We have chosen so that 



3(2 -7)M f< /3 + 2cj(0 



?0 



(43) 



The integration constants ai thus characterize the initial shear and obey 

E-. = 0, E-.— (44) 

i i 

We can use these results to re- write the equation of motion for the Brans-Dicke 
field, Eq. (E5), solely in terms of (f) and our time coordinate ^: 



+ 



4-37 



+ 



(4 - 37)e 



1 



(45) 



The solutions to the Bernoulli equation, Eq. (^), can be cast into the particularly 
simple form 

In f ^) = (4 - 37) f 4r\^i , (46) 
where 0o is an integration constant, by absorbing z(^) into another function. 



AO ^4-37 
M2 



-90 



with go another integration constant. 

Comparing the expression for (j)'/(f) obtained from Eq. 
Eq. (|l|) shows that 



(47) 



with that from 
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(G0) 



(37-4)/2 q3(2-7) 



{G(p)a 



3(2-7) 



M 



M 



(4J 



In other words, 4M'^g{C,) is just the energy density term on the right-hand side 
of the constraint Eq. (^OD , and the definition of g in Eq. (|47|) is precisely this 
constraint equation, where 



9o 



4M2 



(49) 



Thus g{^) must be non-negative. 
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Note that the behavior of the individual scale factors in the string frame, 
aiiO: is given by 



'Mi 



3(2-7) 

exp 



2MJ g{^) 



(50) 



which gives an analytic expression dependent on our ability to perform the inte- 
gration in the exponential. Note that we always have the shear oc (G0)"^a~^, 
where the average scale factor obeys 

.-.^g). (51) 

The shear term, the constant Qq, in the constraint Eq. (^0]), or equivalently 
Eq. (^7[), can only remain important if both g and y"^ remain bounded; the latter 
case implying that remains bounded. Whenever ^ — > cx) the model isotropizes. 

Note that singularities in the Jordan frame, by which we mean here points at 
which the volume scale factor, a^, vanishes, can occur only ii g or (p ^ oo. 
In fact the latter case can be shown also to require that g/z^ and thus 
the dynamical role played by matter (excluding stiff fluid with 7 = 2) at the 
singularity is negligible in anisotropic models. Thus we reserve a discussion of 
the nature of the singularity until Sect. (|5.2|) where we discuss vacuum and stiff 



fluid models. 

The behavior of the coupling u;(0) which defines the scalar-tensor theory is 
given in terms of z{^) by 

3 + 2u m)] = [z' f . (52) 

The uj{(f)) dependence is only obtained after we have solved for the evolution of 
and u as functions of ^, if we can invert Eq. (^) to find ^(0). In practice a theory 
is chosen by specifying g{^) as a generating function from which 0(^) follows by 
Eq. (H), and uj{^) from Eqs. (|3) and (|2|). We have a(0 from Eq. (|T]), as well 
as aj(^) from Eqs. (^OD, and we can relate our time coordinate ^ to the proper 
time in the Jordan frame, t{^), from Eq. (0). 

4.1 Brans-Dicke gravity 

Brans-Dicke theory is recovered when a; = cuq is a constant and we see from 
Eq. (|52|) that this implies that 

(53) 

In the isotropic case where the shear = 0, setting ^0 = corresponds to what 
is often referred to as the solutions being matter (rather than (p) dominated at 
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early times. In fact because both g and then evolve as the relative terms 
on the right-hand side of the constraint Eq. (^) are strictly proportional and so 
this could more accurately be described as a scaling solution. Clearly all Bianchi 
type I solutions approach this behavior, with a oc and oc ^'^ ^ where 

^ 2(3 + 2a;)(2-7)-2(4-37) 
^ 3(3 + 2cu)(2- 7)2- (4-37)2 ' ^ ' 



4(4 - 37) 



3(3 + 2cj)(2- 7)2- (4-37)2 



(55) 



as ^ ^00, giving Nariai's isotropic power-law solutions, a oc t"^ and cx t", 
where 

m = 2(3 + 2^)(2-7)-2(4-37) 

3(3 + 2cj)7(2 - 7) - (4 - 37) (37 - 2) ' ^ ' 

n = 4(4 - 37) 

3(3 + 2cj)7(2 - 7) - (4 - 37) (37 - 2) ' ^ ^ 



upon integration of the time transformation, Eq. (|5^, irrespective of the initial 
shear. 

However, the ^ — > 00 limit is never reached (for positive 0) for u < 2(7 — 
5/3)/(2 — 7)2 as in this case g ^ ai finite ^ (and would become negative as 
^ 00). Nonetheless this universe does correspond to an infinite proper lifetime 
in the Jordan frame, but with a late time expansion driven by both shear and 
the Brans-Dicke field density, i.e., the model does not isotropize. 



4.2 Approach to general relativity 

Among the wider class of scalar-tensor gravity theories Brans-Dicke behavior 
looks atypical. It only occurs at late times where g{^) oc ^'^ in the limit ^ ^ 00. 
Otherwise the generalized Friedmann constraint equation becomes scalar field 
dominated if g{^) < a^2 for any constant a as — *• cxo, or matter dominated if 
g{^) > a^2 as ^ ^ 00. In this latter case we see from Eq. (^) that the Brans- 
Dicke parameter u must diverge, and [from Eq. (^6])] the Brans-Dicke field (p tends 
to a constant value, in other words we recover the general relativistic behavior at 
late times. 

Recently Damour and co-authors [^, ^ have argued that the general rela- 



tivistic limit acts as a cosmological attractor within the parameter space of more 
general scalar-tensor gravity theories. This occurs when the Brans-Dicke parame- 
ter diverges. In the notation favored by Damour and co-authors this corresponds 
to the Brans-Dicke field (p = F{(p) as defined in Eq. (P having a local maximum 
with respect the field ip (with 7^ 0). To see how this emerges in our notation, 
and in the Bianchi type I cosmologies, we will consider the limit 

F{^) = 00 - Ik^' + 0{^') . (58) 
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This simply corresponds to a pole in the function uj{(j)) 



2u;(0) + 3 



1 

2k 



1 + 



(59) 



This type of behavior occurs as ^ ^ oo and we consider a generating function 

9(0 = gne^" + oie-"-') , (60) 



with n > 1. Equation (p6D then gives the evolution of the Brans-Dicke field as 

4-37 



^ 2(n-l)<7„e{-i)+°^^ 
which, via Eq. (|52|) , gives the above uj{(j)) behavior for 

3(2-7)2n- 1 



-2r^.+l^ 



(61) 



k 



4(4 - 37) n2 



(62) 



Note that there is actually a upper bound on k < 3(2 — 7)^/16(4 — 37). This 
corresponds to the condition that the Brans-Dicke field is over-damped and ap- 
proaches 00 monotonically. For larger k it will execute damped oscillations about 
00 0. 

Using Eq. ( |5TD to find the limiting behavior of the average scale factor and 
using this to determine the relationship between the time coordinate ^ and the 
proper time in the Jordan frame, Eq. (P^ ) , we find (for 7 7^ 0) 



a oc ^2/37^ 
^ ^ ^-2(n-l)(2^7)/n7 



OC a 



-3(n-l)(2-7)/n 



(63) 
(64) 

(65) 



For 7 = both ^ and thus a grow exponentially with respect to the proper time 
t leading to de Sitter expansion as ,^ cxd. 

Note that the shear must vanish relative to the all the other terms in the 
constraint Eq. ( ^7D as we approach general relativity. As ^ ^ 00 we find p oc 
but oc t-^/''. 



5 Massless fields in anisotropic cosmologies 

We will in this section restrict ourselves only to vacuum or to matter consisting 
of the short-wavelength modes (radiation, 7 = 4/3) and long- wavelength modes 
(stiff fluid, 7 = 2) of massless fields |16[. Here we can integrate the equations of 
motion in the Einstein frame without specifying a; (</>). 

For radiation the fluid is conformally invariant and in this case, as we have 
seen, a perfect fluid, i.e. non-interacting, in the Jordan frame remains a perfect 
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fluid in the Einstein frame. On tlie otlier liand a non-interacting stiff fluid in the 
Jordan frame does not remain a perfect fluid in the Einstein frame, but we can 
still deal with the dynamics in this case as the homogeneous Brans-Dicke field 
also becomes a stiff fluid, ip, and although there is an interaction between the 
two components their combined dynamical effect is the same as that of a single 
perfect stiff fluid. Thus if the stiff fluid in the Jordan frame is a homogeneous 
minimally coupled scalar field (/?, we can define the composite scalar, x, by 

dx^ = d^"^ + G(j)dip'^ , (66) 

which obeys the equation of motion for a homogeneous minimally coupled field 

^ + ~e& = 0. (67) 

dt^ dt ^ ' 

The corresponding energy density in the Einstein frame is 

where A is a constant of integration. In the absence of a second field (/? then we 
simply have x = '4'- 

Thus with the conformally transformed matter density and pressure, for ra- 
diation and or stiff fluid in the Einstein frame 

3 ^ -1/4/3 -1/2 ' ^ > 

SttG ^ r MGd) , , 

3 3\/4/3 Y2 ' ^ ^ 

where F and M are non-negative constants, we can obtain the energy density of 
the Brans-Dicke field in the Einstein frame as 

. . 1 fdipV 3 ^2 -4MG0 

P = '=2[-dl) =8^ • ^''^ 

Clearly we require 4MG0 < A"^ for this to correspond to a non-negative energy 
density. Thus in the presence of a stiff fluid, M ^ 0, this places an upper bound 
on the value of the Brans-Dicke field, or equivalently a lower limit on the effective 
gravitational coupling constant in the Jordan frame, G^s ~ > {4M/A'^)G. 

The evolution of the volume scale V{t) can be given if we can solve the Einstein 
equations (|20|) and ( pi]) with total energy density 

SttG A^ T 

and pressure 
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as well as the shear and spatial curvature which will depend on the Bianchi class. 
This Einstein frame solution for radiation and a stiff fluid, in a particular Bianchi 
or Kantowski-Sachs metric, is independent of uj{(f)) and describes the Einstein 
frame behavior for any scalar-tensor theory. 

It is only at the final stage that we must specify uj{(f)) in order to invert 



and thus obtain the evolution of (f){x), which is both the Brans-Dicke field and 
the conformal factor relating the Einstein frame to the Jordan metric. We see 
that the effect of a stiff fluid in the Jordan frame (i.e., M 7^ for a given value 
of A^) is to alter the relation between x and (p. But this is the role, in vacuum, 
solely of and so any scalar-tensor theory, defined by Lj{(f)), plus a stiff fluid in 
the Jordan frame is equivalent to the an effective theory with ujva,c{4') in vacuum 
given by 

3 + 2^.,,(0) ^ ^2 SmG4> + '^^^'^^^ ■ ^^^^ 
However this equivalence is broken by the presence of any matter other than 
radiation (i.e., any matter with non-zero trace of the energy-momentum tensor) 
which interacts in the Einstein frame with the field ip, rather than x- 

Even without knowing the specific form of a;(0) we can note a few general 
features of x(0)- In particular if x diverges this must correspond to either a; ^ 00 
or ln(G0) — s> ±00 (with only 0^0 possible for M 7^ 0). The former case may 
occur as approaches a finite, non-zero constant value in which case the the 
Jordan frame coincides (up to an arbitrary constant factor) with the Einstein 
frame and we must have general relativistic behavior. The latter case corresponds 
to singular behavior of the conformal factor so we may expect radically different 
behavior in the Jordan frame from that in the Einstein frame. On the other hand 
as X — >const we must have either uj —>■ —3/2 or (p ^constant. Thus for any 
u > —3/2 we must also recover general relativistic behavior in this limit. 

5.1 Scalar-tensor theories 

We will give here only three particular examples for which we can perform the 
inversion to find 4>{x) analytically in Eq. (|7^, though we can solve for any uj{(f)) 
using numerical integration. 

5.1.1 Brans-Dicke gravity 

When uj = uq = const ant we find 



1^X(0) = ±J^^lnG0 forM = (76) 
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\A\ + ^/A' - AMG(j)' 
\A\ - y/W^AMG(l) 



forM^O. (77) 



Notice that while the presence of a stiff fluid leaves the evolution of the scale 
factor in the Einstein frame unaltered (i.e. independent of the value of M for a 
given A'^) it affects the form of x(0) and thus the evolution of (p and the conformal 
transform back to the Jordan frame. 
Inverting these relations gives 



G0(x) = exp{Po\H^x] forM = (7J 



3 
A^ 



4Mcosh2 (^f^Jl^x 

as shown in Fig.(l), where 



for M ^ (79) 



Thus at an initial singularity in the Einstein frame where x oo we must have 
— > or ^ oo for M = 0. In the presence of a stiff fluid (M ^ 0) only — 
is possible. 

5.1.2 uj{(p) with Brans-Dicke and G.R. limits 

A alternative choice of uj{(j)) that displays both a Brans-Dicke regime and general 
relativistic behavior is 

2^(0) + 3 = (2^0 + 3)7-^^ . (81) 

(0 - 0O)^ 

Clearly for > 0o we have uj > ujq and uj approaches this lower limit as 00, 
while as — > 00 we find uj 00. 
Equation ( |74D then yields 



. J?^,Jt^] faM^O, (82, 



l2uJo + 3 A^ 



3 V - 4MG0 

Note that x ^ 00 as — *• 0o, and also as — >• 00 when M = 0. But because 
< A'^/AMG in the presence of the stiff fluid the latter Brans-Dicke limit cannot 
be reached. 
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This can be inverted to give 

= 00 



1 + exp I \I^^PoX 



for M = 



ImG ~ V 4MG ^ 
where /5o was defined in Eq. (^0]) and 



tanh- I -Bo]l X 



(84) 

for M ^ (85) 



Bo = Po\ 



'A^ -AMG(f)o 



A^ 



(86) 



The function 0(x) is plotted in Fig.(l) 



5.1.3 Axion-dilaton string cosmologies 

While the effective action derived in the low energy limit of string theory is 
sometimes referred to as Brans-Dicke gravity with u = —1, this is strictly only 
true when all other matter fields are minimally coupled to the Jordan, or 'string', 
metric. This is generally not true in string theory, but we will show that at least 
for other massless fields such as the antisymmetric tensor field, H^^x, and moduli 
fields, appearing in the low energy effective action, the techniques developed 
for scalar-tensor gravity can also be applied to string theory. (For a more detailed 
discussion in the case of isotropic FRW cosmologies see p3| .) 

The background field equations of motion reduced to four dimensions can be 
derived from the action W\ 



S 



167rG 



-ge 



1 



12 



^7) 



Conventionally the variable gravitational coupling is represented by the dimen- 
sionless dilaton field ip, which is simply related to the Brans-Dicke field = 
l/(G'e'^). If the antisjTiimetric tensor field is a function only of the 4-dimensional 
spacetime coordinates then it can be described by a single (pseudo-) scalar "axion" 
field h, where 



,K 1 



and e^j^AK is the antisymmetric volume form. 

Conformally transforming to the Einstein frame the above action then be- 
comes 



S 



1 



59) 



Thus we see that the dilaton is simply related to our field defined in Eq. (|TT]), 
as dijj = —dcp/VlQirG. For homogeneous fields we can also define the composite 
scalar field xiv) such that 



dx = # + 



n 



d(3' + 



levrG 



-Ah' 



(90) 
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whose total energy momentum tensor is conserved. Although the axion field in 
particular is coupled directly to the dilaton, we can still integrate the equations 
of motion to deduce the relative energy densities in the Einstein frame: 



Px = 
Ph = 



1 (dx\ _ 3 

2 \di) ~ 8^4K2 ' 

n fd^ _ 3 B '^ 

1 (dh\ 



327rG \ dt 



SnG t/2 



and thus 



1 /#\ _ 



3 - - AMe-^'t' 



SttG 



y2 



(91) 
(92) 
(93) 

(94) 



Once again we see that the presence of other massless fields places an upper 
bound on the value of the Brans-Dicke field, or equivalently a lower bound on the 
dilaton, 



-<p < 



- 52 



AM 



The relation between x and (p is then 
' A 



'leTrG 



2VSJ VA^ - 52 - 4Me-2^ 
1 



d(p , 



A2 



1 I I2 



for M = , 



(95) 

(96) 
(97) 



2V3V^'-52 

^ ^ V V^2 _ 52 _ ^^2 _ 52 _ 4^6-2^=^ 



for M =^ 0(98) 



Inverting this gives 



exp I 2V3i/^^— — ^W^^^x 



A2 



2M 



>12_52 



cosh 



/A2 - ^2 



A2 




\/647rG'x 



For M = and 5 = this coincides with the Brans-Dicke result with u = —1. 
For M — but B the effective Brans-Dicke parameter lies in the range 
— 3/2<a;<— 1. In the presence of the antisymmetric tensor field (M 7^ 0) the 
form of = {Ge'^)~^ differs from a purely Brans-Dicke result. 
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Thus we can use the same results as we will use for general scalar-tensor 
gravity models with stiff fluid (with or without radiation) in the Einstein frame 
to derive the general solutions to the low energy string effective action including 
homogeneous antisymmetric tensor and moduli fields (with or without radiation). 



5.2 Anisotropic models 
5.2.1 Bianchi type I 

Considering again the spatially fiat, anisotropic metric given in Eq. (p5D, we will 
introduce the variable X = = {G(f))a'^ and the conformally invariant time 
coordinate drj = dt/a = dt/a. Remembering that the shear in the Einstein frame 
is given by cf^ = 32^/45^, and using p as given in Eq. (|7lD , so that 



3 i^^^i 



V 3 + 2cj X 

and p given in Eq . (|69|) , we find the constraint Eq. (|20|) , in vacuum or with a stiff 
fluid in the Jordan frame (but for F = 0), can be simply written as 

X'2 = ^2 + e2 . (102) 

This is precisely the same constraint equation as solved in the case of scalar-tensor 
gravity in fiat FRW models in vacuum [15| or with stiff fluid [16 and we have 



X = a^ = VA^Tt? \7] - r]o\ . (103) 

In the presence also of radiation in the Jordan frame the constraint Eq. (pO]) 
becomes 

= ^2 + ^2 + 4rX , (104) 

and so 

X = d^ = \ri - rio\ (VA^TT^ + T |r/ - r/o|) . (105) 



These are well-known general relativistic results for the evolution of Bianchi 
type I models in the presence of radiation and stiff fluid. As one would expect 
the stiff fluid and shear dominate the evolution near the singularity (a 0) but 
the radiation term will dominate as a — > oo. 

The behavior of this averaged scale factor, together with the shear 

5' = , , (106) 

4|r/-r/o|3(VA2 + s2 + r|r/-r/o|)3 

describe the general evolution. However the Bianchi type I metric has three 
degrees of freedom so there is a degeneracy within the evolution we have described 
so far depending on how the expansion and shear is divided between the three 
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scale factors. Solving the equations of motion for each scale factor in the Einstein 
frame we find 

a^ = \v-Vo\"''^'{VA^TT^+T\r]-Vo\y'^"''^'^ , (107) 

where the definitions of the overall expansion and shear give two constraints on 
the three new integration constants: 

Thus although there is a unique (isotropic) late-time behavior in the the presence 
of radiation (F 7^ 0), where Sj oc \ri — ?7o|, the initial behavior (for |?7 — r/ol 
\/ A"^ + S^/r) is dependent upon the choice of integration constants, Cj. [See 
Fig. (2)]. 

This early power-law evolution of the scale-factors includes the pure general 
relativistic "Kasner" vacuum solution with A = 0. Although the averaged scale 
factor, a, always vanishes at t]q, this singularity can be point-like (all Cj > for 
all i) or linear (only one Cj < 0) unlike the Kasner behavior for which only the 
linear singularity is possible (except when Cj = 1 and Ci-^j = and the singularity 
is planar). Indeed for A"^ > only the point-like singularity is possible [^. 

To recover the full scalar-tensor results in the original Jordan frame we must 
calculate x as a function of time from Eq. 



IQttG A 

-^ixiv)-Xo) = ^^^^^^ ln|^-^o| forr = (109) 

A 1^ r\v -V o\ forr^OllO) 



VW+T^ VW+T^ + T\r] - r]o\ 

and then use 4>{x)^ dependent on the form of uj{(f)), to specify the evolution of the 
Brans-Dicke field, and thus the conformal factor relating these solutions back to 
the Jordan frame variables. Note that the field x must always diverge as — 779. 
(In the absence of radiation it also diverges at 77 — >• ±00.) 

As noted earlier, the definition of x(0) Eq.([7^) requires either u; — 00 or 
ln(G0) ±00. The former case where tends to a finite, non-zero value of will 
lead to a purely general relativistic result as X — > 0. In the latter case we can 
expand about the point X —>■ and we find that in the Jordan frame we have 

<p^\v-vor, (111) 

a oc |r/-r7o|('~"^/' , (112) 

a^oclv-Vof''""^^' , (113) 
t(x\r]-7]o\^^-"^/\ (114) 



if — s> uj{0) 7^ —3/2 and n 7^ 3, where from Eq.( |101| ) 
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Thus without specifying the full 0^(0) we can deduce a number of features of 
the possible behavior in the Jordan frame: 



1. In the presence of shear and/or a stiff fluid in the Jordan frame, a "bounce" 
(where the volume factor V is stationary, dV/dt = 0) will occur whenever 

^=if£:-£Uo. (116) 

which, from Eqs. ( [101| ) and ( |104| ), requires 



3 , 3 f A'-4MG<p \ 



We thus see that a bounce is only possible when 

3 



2<^<0. (118) 



In a fiat FRW metric (S^ = 0) we require 00 = [T^. The presence of shear 
or a stiff fluid requires a negative value of u, producing a negative effective 
energy density for the Brans-Dicke field in the Jordan frame, to produce a 
bounce. For u;(0) less than that given in Eg. (|117|) , the averaged scale factor 



in the Jordan frame actually grows as both X and approach zero. 

2. The singularity at X = is always present in the Einstein frame, and is 
bound to be anisotropic for 7^ 0, with — 00 as X — >■ 0. However 
we noted in Sect. (0) that in the Jordan frame the Brans-Dicke field exerts 
an anisotropic pressure proportional to the shear and one might wonder 
whether it is possible for this to suppress the shear as one approaches the 
singularity. This is indeed possible. The shear in the Jordan frame in a 
Bianchi type I model is given by 

2 3E2G0 , , 

so if vanishes faster than the anisotropic initial singularity in the 
Einstein frame becomes isotropic in the Jordan frame. Again from Eqs. ( |7^ 
and (|104| ) we find that this requires — > and 



3. Is it then possible that the singularity as X — »■ in the Einstein frame can 
be avoided completely in the Jordan frame? First note that the condition 
that da/dt = [Eq. ( |117| )] is incompatible with the condition that the shear 
should vanish [Eq. (|12CI|) ], or even remain finite, as X and both approach 
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zero. However, the expansion 9 = 3{da/dt)/a can remain finite even though 
da/dt diverges if the average scale factor a grows fast enough. Using the 
above results for a general scalar-tensor theory in a Bianchi type I metric 
with a stiff fluid, we can write the expansion as 





l±x-^ ^\ 7-. ^ . fl211 



g2 _ 3(^2 + S 
4 

We see that whenever (j) a.s X ^ such that the shear o"^ ^ 0, then 
the expansion will also vanish (for uj ^ —?)/2). This also coincides with the 
case where f] ^ r]o takes an infinite proper time in the Jordan frame. Thus 
although these models reach an anisotropic singularity in a finite time in the 
Einstein frame, this corresponds to a non- singular, shear-free infinite proper 
lifetime in the Jordan frame. The condition for this to occur is simply the 
condition for the shear to vanish as X ^ given above, which can occur 
for 

3 4 
-- <uj < — , (122) 
2 3' ^ ' 

m a fiat FRW metric (S^ = 0), or for the more limited range given in 
Eq. ( |120[ ) depending on the relative strength of the anisotropy.0 Remember 
that = coincides with x(0) ~^ cxo for uj ^ —3/2 and so — can be an 
attractor solution as X ^ where we have shown x must diverge. 

4. If grows sufficiently rapidly as we approach X ^ then the dynamical 
effect of other matter, neglected here, may no longer remain negligible. 
The condition for the density of isotropic matter with barotropic index 7 
to decrease with respect to the shear or Brans-Dicke field density in the 
Einstein frame as r] ^ rjo is that 

For smaller values of uj it may no longer be possible to neglect the effect 
of this matter as a — > 0. In particular in the case of the non-singular 
solutions given above with vanishing shear and expansion (in the Jordan 
frame) as rj ^ rjo, the relative density of barotropic matter always grows 
as we approach rjo for fluids with 7 < 1. We have also neglected here any 
possible anisotropy in the matter content. This, of course, may well play 
an important role in anisotropic solutions, but is beyond the scope of this 
paper and we leave this for future work. 



■^The requirement that lu < —4/3 to avoid the singularity was pointed out by Nariai in 
the case of Brans-Dicke gravity. 
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5.2.2 Bianchi type V 

The next simplest case to consider is that of a Bianchi type V metric which can 
be written as 

ds^ = -de + a'' it) (dx^ + e^" [L\t)dy^ + L-^{t)dz^] ) (124) 

Just as Bianchi type I includes the flat FRW model, Bianchi type V includes 
the open FRW model in the isotropic limit, L =const. Here we have already 
introduced the averaged scale factor a, so that the expansion is simply 

a dt 

and the shear is then 

IdL , , 

Like open FRW models, the homogeneous hypersurfaces have a negative spatial 
curvature 

^'^R = ^ , (127) 

but because this is a function of the averaged scale factor, i.e. does not select out 
a particular direction, it does not drive the shear. 

Thus transforming to the Einstein frame where the Brans-Dicke field decouples 
from the spacetime curvature, we find that the shear again evolves as a free field, 

= 3S^/4a^, just as in Bianchi type I. Introducing X = = G0a^ and 
dr] = dt/d = dt/a as before, the constraint equation, including stiff fluid and 
radiation, simply becomes 

X'2 - ^ ^2 ^ ^2 ^ , (128) 

This is mathematically identical to the constraint equation solved in the case 
of scalar-tensor gravity for an open FRW model [jl6l. We can integrate this to 
obtain 



X = d' = ^ , (129) 

1 — 

where we have written 

r = tanh \r] — rjol . (130) 

This variable, r, turns out to play much the same role in the presence of negative 
spatial curvature as the conformal time, rj, does in the spatially fiat case, as is 
the case in FRW models |16]. 

At late times [i] ±oo) r ^ 1 and the evolution becomes curvature domi- 
nated, as we would expect in a non-infiationary universe. Because the curvature 
does not drive the shear, 

O- = : — , (131) 



4r3 (VA2 + S2 + Ft 
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it vanishes as r — * 1. 

At early times r ~ |?7 — ?7o| ^ 1 and so the curvature is irrelevant and we 
recover a Bianchi type I solution. However unlike the general Bianchi type I, the 
metric has only two degrees of freedom and so its evolution in the Einstein frame 
is completely described by the expansion and shear. Integrating the expression 
for a gives 

= Llr^^^/-^^^' forr = 0, (132) 

The behavior near the singularity in Bianchi type V metrics is thus only a subset 
of the Bianchi type I solutions, given in Eqs. ( p.07| ) and ( fLOSp , restricted to C2— C3 = 
\/3E/(A^ + T?). Thus (for a given A and V) they are parametrized solely by the 
choice of S. 

Similarly, because the Brans-Dicke field is decoupled from the spacetime cur- 
vature in the Einstein frame, the field 

A 

Xo) = Inr forr = 0, (134) 

A Vt 

In ^ for r 7^ , (135) 




and thus 0(x) approaches a fixed value 0oo as r 1 (unless x(0) is singular at 
this point). Thus the evolution of the metric in the original Jordan frame will 
also approach that in the Einstein frame at late times. 

5.2.3 Bianchi type III and Kantowski-Sachs 

Here we will write the metric in the Jordan frame as 

ds^ = ~dt^ + al{t)dx^ + al{t)(dy'^ + s^{y)dz^) , (136) 

where 

{siny for Kantowski — Sachs, 
y for LRS Bianchi I, (137) 

sinhy for Bianchi HI. 

The volume scale factor V = aia\. If we introduce a conformally invariant time 
coordinate d^ = ^ = ^ and let X = dia2, we can write the Einstein equations 
as 




X" + 4kX = 87TG{p~p)Xal (138) 
= AnG{p-p)dl (139) 

= 87iGp + pal, (140) 
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where ' = d/d^ and k = +1,0, —1 corresponds to Kantowski-Sachs, Bianchi I or 
Bianchi III respectively, in analogy with FRW models. 

We have only been able to solve these equations analytically in the presence 
of a stiff fluid plus the Brans-Dicke field. Even in the case of radiation, where 
one can obtain both p and p as functions of Si and 02, the resulting equations 
for Si and S2 still cannot be integrated. However, we see that a stiff fluid in the 
Einstein frame, like the Brans-Dicke field, does not enter the first two equations 
and so (for F = 0) we can integrate both of these directly to give 

X = ^s{2[^-^o]) , (141) 

with the appropriate function s{x) defined above and thus 

~ ~ c 

Oi = Cti^<T , 

(143) 

where we have defined 

I tan(^ — ^o)| for Kantowski — Sachs, 

- ^ol for LRS Bianchi I, (144) 

I tanh(^ — ^o)| for Bianchi III. 

A, C and Si,,, are integration constants and S2* = A/ai^. The evolution of the 
two scale factors is shown in Fig. (3). 

The energy density of the stiff fluid in the Einstein frame is given by Eq. 
where the constraint Eq. ( |140| ) requires 



A' = \dlal{l-C^). (145) 

and so clearly we must have < 1 in the presence of a stiff fluid in the Einstein 
frame, or C = ±1 in general relativistic (i.e. (j) =const) vacuum. 

Note that we again obtain Kasner type solutions near the singularity at ,^ ^ 
^0, independent of the sign of k. The curvature becomes irrelevant and we recover 
a LRS subset of the Bianchi type I solutions, a, oc |?7 — //ol^'^'''^! where drj = 
dt/{didiy/^ and 

C 1 — C 

ci = j—^ and C2 = C3 = j—^ > . (146) 

These q obey the relations given in Eq. ( |108| ) with 

^2 f 2aua2.{l-2C) V _ ^^^^^ 
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The area of the two-dimensional sub-space always vanishes as we approach the 
singularity, while the remaining spatial dimension is free to diverge (for C < 0) or 
collapse (C > 0). Thus the singularity is linear or point-like. In the Kantowski- 
Sachs case we find another singularity at ^ ^ + ^/2- This is like the initial 
singularity but with C — — C, and thus a point-like singularity is followed by a 
linear singularity, or vice versa (unless C = in which case Si remains constant 
throughout). 

We can then integrate Eq. (^S|) to obtain 

-X-o) = (148) 




^ /i^.n^. (149) 

Thus X must diverge both at early and late times in the Kantowski-Sachs (or 
Bianchi type I) metric, but the expanding universe becomes curvature dominated 
in the Bianchi type III metric with positive spatial curvature and x coasts to a 
fixed value x ^^4(1 - C2)/3. 

5.2.4 LRS Bianchi type IX 

To give an example of an anisotropic cosmology with closed spatial hypersur- 
faces we consider a Bianchi type IX metric |^ 0, whose homogeneous spatial 
hypersurfaces have volume V = 1671010203, so the expansion 

6 = ^ + ^ + ^. (150) 

Ol 02 Ct3 

When di/oi = 02/02 = 03/03 we recover the closed FRW model. 

As the curvature terms in anisotropic models become more complicated, our 
ability to give analytic solutions becomes more restricted. In the case of Bianchi 
type IX we can only give analytic solutions for the locally rotationally symmetric 
(LRS) case (where 02 = 03) plus stiff fluid in the Einstein frame corre- 



sponding to vacuum or stiff fluid solutions in scalar-tensor gravity. This will not 



show the chaotic behavior of the more general Bianchi IX metric nor the 
isotropizing effect of matter such as isotropic radiation. Such issues would re- 
quire a phase-space analysis and is beyond the scope of this paper but is perhaps 
a topic worthy of a investigation in its own right. 

The equations of motion for the scale factors in the Einstein frame are then 



2a; 




(151) 



2o| ' ^'''^ 
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plus the constraint equation 



'lU2?l| + ^^11^ = 1^. (153) 



where the only energy density is p-^ given in Eq. (p8|). 

Introducing the volume weighted time coordinate = dtjiaia^ and using 
the variables x = Aaf and y = ^ala^, we can re- write the equations of motion as 

(154) 

(155) 

which can immediately be integrated to give 

~2 ^1 




cosh Wi(^ — C- 



1) 

~al = ^2^"^h^i(C-Ci) 
wi cosh W2(C - C2) 

as shown in Fig. (4), where Ci, C2, and W2 are integration constants. Clearly, wi 
must be positive and we can also take W2 to be positive without loss of generality. 
The constraint equation requires 

= Awl ~ ^1 ■ (157) 

This is sufficient to ensure that Wi < 2w2 and both 0,2 and Si approach zero 
as C ^ ioo (for non-zero Wi and W2), the interval between this big bang and 
big crunch taking only a finite proper time in the Einstein frame. Near these 
singularities the scale factors evolve as power-laws, with respect to proper time 
(or conformal time), and we recover another one-parameter (for given A) subset 
of the Bianchi type I solutions, with 

wi 2wi-wi 
ci = , C2 = C3 = . 158 

AW2 — Wi AW2 — Wi 

The initial (and final) shear then approaches 3S^/4a^ with = 4(^2 — 

wif/9. 

Note that the equation for the scalar field energy density, Eq. (|68|) , shows that 
X is just proportional to the volume weighted time, (, and so we have 



\/^(x-Xo) = ±AC, (159) 
which must also diverge as we approach both the initial and final singularity. 
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6 Conclusions 



We have derived a number of new exact solutions for anisotropic cosmologies in 
scalar-tensor gravity theories where previous results were in general restricted 
to isotropic models |1^, |I6| or the particular case of Brans-Dicke gravity 



|T0|, |T^, ^ |2^, |T2|. In the original, or Jordan, frame the non-minimal coupling 
of the Brans-Dicke field introduces an effective anisotropic pressure proportional 
to the existing shear. This can significantly modify the evolution even in vacuum, 
especially near a singularity. It is possible for the Brans-Dicke field to reverse the 
collapse as we approach an anisotropic singularity and lead to a non-singular, 
isotropic expanding universe. In the Brans-Dicke gravity theory the deviation 
from general relativistic behavior is strictly limited due to present day obser- 
vational limits on the Brans-Dicke parameter uj However in more general 
scalar-tensor gravity the present day value need not constrain the value of uo near 
an initial singularity. 

We have emphasized the importance of using the conformally related Einstein 
metric where the Brans-Dicke field is minimally coupled with respect to the metric 
and so the usual general relativistic results hold good. Here, for instance, there is 
no anisotropic pressure due to the Brans-Dicke field. Non-Einstein behavior only 
appears in the transformation back to the Jordan frame. 

The analytic complexity of the evolution in the Einstein frame lies in the in- 
teraction introduced directly between the Brans-Dicke field and ordinary matter. 
Except in the particular cases of vacuum, radiation or a stiff fluid, this forces us 
to restrict our results for general barotropic fluids [p = (7 — l)p] to the spatially 
fiat Bianchi type I metric. Here the solutions for general (^(0) theories are charac- 
terized by a generating function g{C}i where ^ is a time coordinate. Brans-Dicke 
gravity corresponds to the case where oc + ^o)^- 

The shear vanishes relative to the density of matter and/or the Brans-Dicke 
field (i.e., the models isotropize) whenever |^| 00. This requires simply that 
g remains positive, or equivalently that the Brans-Dicke parameter uj > 2(7 — 
5/3)/(2 — 7)^, as 1^1 ^00. If g < a^"^ for any constant a > as ^ — cxd then 
the dynamics become dominated by the Brans-Dicke field density at late times. 
Conversely, if g > a^"^ then the models become dominated by the barotropic 
matter and we recover the general relativistic behavior where the Brans-Dicke 
field (j) ^const and a; — >■ cxo. Brans-Dicke gravity is seen to correspond to the 
particular case where the relative densities of the Brans-Dicke field and ordinary 
matter remain proportional. 

Singularities in the Jordan frame (where the volume scale factor vanishes) 
occur only when the effect of ordinary matter (with p < p) becomes negligible 
compared with the shear and Brans-Dicke field energy density. Thus it is sufficient 
to consider only the vacuum or stiff fluid {p = p) models to discuss the approach 
to the singularity. 

We have shown that in the presence only of radiation and a stiff fluid (equiva- 
lent to the short and long wavelength modes of massless fields [0), or in vacuum. 
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the evolution of the scale factor in the Einstein frame is independent of the form 
of uj{(j)) and corresponds exactly to the standard general relativistic evolution 
with radiation and a stiff fluid. This enables us to give results for arbitrary a;(0) 
theories in Bianchi type I, III, V and LRS type IX and Kantowski-Sachs metrics. 
At singularities in the Einstein frame these all approach a Bianchi type I solution 
as the spatial curvature becomes negligible. At this point the energy density of 
the scalar held in the Einstein frame, %, must diverge. The relation between 
this held and the original Brans-Dicke held depends on the form of but 
in any case a necessary condition for the divergence of x is that ln(0) diverges 
or C(j — s> oo. However the simplest general relativistic limit, considered above as 
a late-time attractor in Bianchi type I, in vacuum or in the presence solely of 
radiation or stiff fluid, is not suflicient to lead to a divergence of x ^-nd does not 
act as an attractor solution at the initial singularity. On the other hand — > 
is an attractor which can lead to non-singular behavior in the Jordan frame. 

Thus while many scalar-tensor gravity theories do approach general relativity 
at late times, where the role of other barotropic matter is dominant and anisotropy 
vanishes, they may nonetheless give markedly non-Einstein behavior near the 
initial big bang where anisotropy is important. 
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Figure Captions 



Figure 1 

The Brans-Dicke field 0, as a function of the minimally coupled field in the 
Einstein frame, x- Brans-Dicke gravity (with uj = 1) is given by solid 
lines (for two choices of sign of (3q) in vacuum or radiation {A^ = 2, M = 0) or 
by dotted line in presence of a stiff fluid (M = 0.2). For co'(0) given in Eq. (^Tj), 
with 00 = 0.8, the dashed line gives the 0(x) for M = 0, or dotted-dashed line 
for M = 0.2. 

Figure 2 

Evolution of the three independent scale factors Sj, in the Einstein frame, in 
Bianchi type I metric as given in Eq. ( |107| ) with = 2 and with radiation 
(solid line, F = 2) or without (dashed line, F = 0). We have chosen ci = —1/4, 
C2 = 1/2, and C3 = 3/4 in both cases. 

Figure 3 

Evolution with respect to time coordinate ^ of the two independent scale factors Si 
and 0,2, in the Einstein frame, in Kantowski-Sachs (solid line) or Bianchi type III 
metric (dashed) as given in Eq. ( |143| ). We have taken C = 0.3. 

Figure 4 

Evolution with respect to time coordinate ( of two independent scale factors Si 
and 0,2, in Einstein frame, in LRS Bianchi type IX metric as given in Eq. (|156D . 
We have taken = 2 and Wi = 5. 



30 



-5 5 

(87VG)^/2 X 



DWANDS 22-Dec-1994 10: 




DWANDS 22-Dec-1994 11:30 




DWANDS 22-Dec-1994 13:52 




DWANDS 22-Dec-1994 14: 8 



